/ Celestin Apprentice 4 / Apprentice-Release4.iso / Languages / RLaB 1.18c / manual / eign.tex

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 73 65 63 74 69 6f 6e | 7b 65 69 67 6e 20 2d 2d |\section|{eign --|
|00000010| 2d 20 4e 6f 6e 2d 73 79 | 6d 6d 65 74 72 69 63 20 |- Non-sy|mmetric |
|00000020| 45 69 67 65 6e 20 44 65 | 63 6f 6d 70 6f 73 69 74 |Eigen De|composit|
|00000030| 69 6f 6e 7d 0d 5c 6c 61 | 62 65 6c 7b 66 6e 3a 65 |ion}.\la|bel{fn:e|
|00000040| 69 67 6e 7d 0d 5c 69 6e | 64 65 78 7b 65 69 67 6e |ign}.\in|dex{eign|
|00000050| 28 29 7d 0d 5c 69 6e 64 | 65 78 7b 65 69 67 65 6e |()}.\ind|ex{eigen|
|00000060| 20 64 65 63 6f 6d 70 6f | 73 69 74 69 6f 6e 7d 0d | decompo|sition}.|
|00000070| 5c 62 65 67 69 6e 7b 72 | 61 69 6c 7d 0d 45 69 67 |\begin{r|ail}.Eig|
|00000080| 65 6e 4e 20 3a 20 27 65 | 69 67 6e 28 27 20 28 20 |enN : 'e|ign(' ( |
|00000090| 4d 61 74 72 69 78 20 7c | 20 4d 61 74 72 69 78 20 |Matrix || Matrix |
|000000a0| 27 2c 27 20 4d 61 74 72 | 69 78 20 29 20 27 29 27 |',' Matr|ix ) ')'|
|000000b0| 20 3b 0d 5c 65 6e 64 7b | 72 61 69 6c 7d 0d 0d 5c | ;.\end{|rail}..\|
|000000c0| 76 65 72 62 2b 65 69 67 | 6e 2b 20 77 69 74 68 20 |verb+eig|n+ with |
|000000d0| 61 20 73 69 6e 67 6c 65 | 20 73 71 75 61 72 65 20 |a single| square |
|000000e0| 61 72 67 75 6d 65 6e 74 | 20 6d 61 74 72 69 78 20 |argument| matrix |
|000000f0| 63 6f 6d 70 75 74 65 73 | 20 74 68 61 74 20 6d 61 |computes| that ma|
|00000100| 74 72 69 78 27 73 0d 65 | 69 67 65 6e 76 65 63 74 |trix's.e|igenvect|
|00000110| 6f 72 73 20 61 6e 64 20 | 65 69 67 65 6e 76 61 6c |ors and |eigenval|
|00000120| 75 65 73 2e 20 54 68 65 | 20 72 65 73 75 6c 74 73 |ues. The| results|
|00000130| 20 61 72 65 20 72 65 74 | 75 72 6e 65 64 20 61 73 | are ret|urned as|
|00000140| 20 61 20 6c 69 73 74 20 | 77 69 74 68 0d 65 6c 65 | a list |with.ele|
|00000150| 6d 65 6e 74 73 20 5c 76 | 65 72 62 2b 6c 76 65 63 |ments \v|erb+lvec|
|00000160| 2b 2c 20 5c 76 65 72 62 | 2b 76 61 6c 2b 20 61 6e |+, \verb|+val+ an|
|00000170| 64 20 5c 76 65 72 62 2b | 72 76 65 63 2b 20 77 68 |d \verb+|rvec+ wh|
|00000180| 69 63 68 20 61 72 65 20 | 74 68 65 0d 6c 65 66 74 |ich are |the.left|
|00000190| 20 65 69 67 65 6e 76 65 | 63 74 6f 72 73 2c 20 65 | eigenve|ctors, e|
|000001a0| 69 67 65 6e 76 61 6c 75 | 65 73 20 61 6e 64 0d 72 |igenvalu|es and.r|
|000001b0| 69 67 68 74 20 65 69 67 | 65 6e 76 65 63 74 6f 72 |ight eig|envector|
|000001c0| 73 20 72 65 73 70 65 63 | 74 69 76 65 6c 79 2e 20 |s respec|tively. |
|000001d0| 49 6e 20 67 65 6e 65 72 | 61 6c 2c 20 74 68 65 73 |In gener|al, thes|
|000001e0| 65 20 77 69 6c 6c 20 62 | 65 0d 63 6f 6d 70 6c 65 |e will b|e.comple|
|000001f0| 78 20 71 75 61 6e 74 69 | 74 69 65 73 2e 20 20 0d |x quanti|ties.  .|
|00000200| 0d 5c 76 65 72 62 2b 65 | 69 67 6e 2b 20 77 69 74 |.\verb+e|ign+ wit|
|00000210| 68 20 74 77 6f 20 61 72 | 67 75 6d 65 6e 74 73 20 |h two ar|guments |
|00000220| 63 6f 6d 70 75 74 65 73 | 20 74 68 65 20 65 69 67 |computes| the eig|
|00000230| 65 6e 76 65 63 74 6f 72 | 73 20 61 6e 64 20 0d 65 |envector|s and .e|
|00000240| 69 67 65 6e 76 61 6c 75 | 65 73 20 6f 66 20 74 68 |igenvalu|es of th|
|00000250| 65 20 72 65 61 6c 20 67 | 65 6e 65 72 61 6c 69 7a |e real g|eneraliz|
|00000260| 65 64 20 73 79 6d 6d 65 | 74 72 69 63 20 28 20 6f |ed symme|tric ( o|
|00000270| 72 20 63 6f 6d 70 6c 65 | 78 20 67 65 6e 65 72 61 |r comple|x genera|
|00000280| 6c 69 7a 65 64 0d 48 65 | 72 6d 69 74 69 61 6e 29 |lized.He|rmitian)|
|00000290| 20 64 65 66 69 6e 69 74 | 65 20 65 69 67 65 6e 70 | definit|e eigenp|
|000002a0| 72 6f 62 6c 65 6d 2e 20 | 54 68 65 20 65 69 67 65 |roblem. |The eige|
|000002b0| 6e 76 61 6c 75 65 73 20 | 61 6e 64 20 65 69 67 65 |nvalues |and eige|
|000002c0| 6e 76 65 63 74 6f 72 73 | 20 61 72 65 0d 72 65 74 |nvectors| are.ret|
|000002d0| 75 72 6e 65 64 20 69 6e | 20 61 20 6c 69 73 74 20 |urned in| a list |
|000002e0| 61 73 20 66 6f 72 20 74 | 68 65 20 66 69 72 73 74 |as for t|he first|
|000002f0| 20 66 6f 72 6d 2e 0d 0d | 54 68 69 73 20 72 6f 75 | form...|This rou|
|00000300| 74 69 6e 65 20 66 6f 72 | 63 65 73 20 75 73 65 20 |tine for|ces use |
|00000310| 6f 66 20 74 68 65 20 6e | 6f 6e 2d 73 79 6d 6d 65 |of the n|on-symme|
|00000320| 74 72 69 63 20 65 69 67 | 65 6e 73 6f 6c 76 65 72 |tric eig|ensolver|
|00000330| 20 2d 2d 2d 20 6e 6f 20 | 63 68 65 63 6b 69 6e 67 | --- no |checking|
|00000340| 0d 69 73 20 64 6f 6e 65 | 2e 0d 0d 0d 0d 5c 70 61 |.is done|.....\pa|
|00000350| 72 61 67 72 61 70 68 7b | 45 78 61 6d 70 6c 65 7d |ragraph{|Example}|
|00000360| 0d 5c 69 6e 64 65 78 7b | 65 78 61 6d 70 6c 65 20 |.\index{|example |
|00000370| 6f 66 20 65 69 67 6e 28 | 29 7d 0d 5c 69 6e 64 65 |of eign(|)}.\inde|
|00000380| 78 7b 65 69 67 6e 28 29 | 21 65 78 61 6d 70 6c 65 |x{eign()|!example|
|00000390| 7d 0d 5c 62 65 67 69 6e | 7b 76 65 72 62 61 74 69 |}.\begin|{verbati|
|000003a0| 6d 7d 0d 3e 20 44 20 3d | 20 72 61 6e 64 28 34 2c |m}.> D =| rand(4,|
|000003b0| 34 29 0d 20 44 20 3d 0d | 20 20 20 20 30 2e 35 31 |4). D =.|    0.51|
|000003c0| 32 20 20 20 20 20 20 30 | 2e 35 32 32 20 20 20 20 |2      0|.522    |
|000003d0| 20 20 30 2e 31 31 31 20 | 20 20 20 20 20 30 2e 36 |  0.111 |     0.6|
|000003e0| 39 33 20 20 0d 20 20 20 | 20 30 2e 38 30 37 20 20 |93  .   | 0.807  |
|000003f0| 20 20 20 20 30 2e 32 35 | 33 20 20 20 20 20 20 30 |    0.25|3      0|
|00000400| 2e 32 31 36 20 20 20 20 | 20 20 30 2e 37 36 32 20 |.216    |  0.762 |
|00000410| 20 0d 20 20 20 20 30 2e | 33 39 37 20 20 20 20 20 | .    0.|397     |
|00000420| 20 30 2e 33 33 39 20 20 | 20 20 20 30 2e 30 32 32 | 0.339  |   0.022|
|00000430| 32 20 20 20 20 20 20 30 | 2e 37 31 37 20 20 0d 20 |2      0|.717  . |
|00000440| 20 20 20 30 2e 33 36 31 | 20 20 20 20 20 30 2e 30 |   0.361|     0.0|
|00000450| 32 35 31 20 20 20 20 20 | 20 30 2e 33 35 36 20 20 |251     | 0.356  |
|00000460| 20 20 20 30 2e 30 36 30 | 36 20 20 0d 3e 20 69 73 |   0.060|6  .> is|
|00000470| 73 79 6d 6d 28 44 29 0d | 20 20 20 20 20 20 20 20 |symm(D).|        |
|00000480| 30 0d 3e 20 65 69 67 65 | 6e 4e 20 3d 20 65 69 67 |0.> eige|nN = eig|
|00000490| 6e 28 44 29 0d 20 65 69 | 67 65 6e 4e 20 3d 0d 20 |n(D). ei|genN =. |
|000004a0| 20 20 6c 76 65 63 20 20 | 20 20 20 20 20 20 20 72 |  lvec  |       r|
|000004b0| 76 65 63 20 20 20 20 20 | 20 20 20 20 76 61 6c 20 |vec     |    val |
|000004c0| 0d 3e 20 44 20 2a 20 65 | 69 67 65 6e 4e 2e 72 76 |.> D * e|igenN.rv|
|000004d0| 65 63 5b 3b 31 5d 20 2d | 20 65 69 67 65 6e 4e 2e |ec[;1] -| eigenN.|
|000004e0| 76 61 6c 5b 31 5d 2a 20 | 65 69 67 65 6e 4e 2e 72 |val[1]* |eigenN.r|
|000004f0| 76 65 63 5b 3b 31 5d 0d | 20 20 20 20 20 20 20 20 |vec[;1].|        |
|00000500| 20 20 20 33 2e 33 33 65 | 2d 31 36 20 2b 20 30 69 |   3.33e|-16 + 0i|
|00000510| 0d 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |.       |        |
|00000520| 20 20 20 30 20 2b 20 30 | 69 0d 20 20 20 20 20 20 |   0 + 0|i.      |
|00000530| 20 20 20 20 20 20 20 20 | 20 20 20 20 30 20 2b 20 |        |    0 + |
|00000540| 30 69 0d 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |0i.     |        |
|00000550| 20 20 20 20 20 30 20 2b | 20 30 69 0d 3e 20 44 20 |     0 +| 0i.> D |
|00000560| 2a 20 65 69 67 65 6e 4e | 2e 72 76 65 63 5b 3b 32 |* eigenN|.rvec[;2|
|00000570| 5d 20 2d 20 65 69 67 65 | 6e 4e 2e 76 61 6c 5b 32 |] - eige|nN.val[2|
|00000580| 5d 2a 20 65 69 67 65 6e | 4e 2e 72 76 65 63 5b 3b |]* eigen|N.rvec[;|
|00000590| 32 5d 0d 20 20 20 20 20 | 20 20 20 20 20 2d 32 2e |2].     |     -2.|
|000005a0| 32 39 65 2d 31 36 20 2b | 20 30 69 0d 20 20 20 20 |29e-16 +| 0i.    |
|000005b0| 20 20 20 20 20 20 2d 31 | 2e 33 39 65 2d 31 36 20 |      -1|.39e-16 |
|000005c0| 2b 20 30 69 0d 20 20 20 | 20 20 20 20 20 20 20 2d |+ 0i.   |       -|
|000005d0| 31 2e 31 31 65 2d 31 36 | 20 2b 20 30 69 0d 20 20 |1.11e-16| + 0i.  |
|000005e0| 20 20 20 20 20 20 20 20 | 2d 39 2e 37 31 65 2d 31 |        |-9.71e-1|
|000005f0| 37 20 2b 20 30 69 0d 3e | 20 44 20 2a 20 65 69 67 |7 + 0i.>| D * eig|
|00000600| 65 6e 4e 2e 72 76 65 63 | 5b 3b 33 5d 20 2d 20 65 |enN.rvec|[;3] - e|
|00000610| 69 67 65 6e 4e 2e 76 61 | 6c 5b 33 5d 2a 20 65 69 |igenN.va|l[3]* ei|
|00000620| 67 65 6e 4e 2e 72 76 65 | 63 5b 3b 33 5d 0d 20 20 |genN.rve|c[;3].  |
|00000630| 20 20 20 20 20 20 20 20 | 20 31 2e 39 34 65 2d 31 |        | 1.94e-1|
|00000640| 36 20 2b 20 30 69 0d 20 | 20 20 20 20 20 20 20 20 |6 + 0i. |        |
|00000650| 20 20 36 2e 39 34 65 2d | 31 37 20 2b 20 30 69 0d |  6.94e-|17 + 0i.|
|00000660| 20 20 20 20 20 20 20 20 | 20 20 20 35 2e 35 35 65 |        |   5.55e|
|00000670| 2d 31 37 20 2b 20 30 69 | 0d 20 20 20 20 20 20 20 |-17 + 0i|.       |
|00000680| 20 20 20 2d 35 2e 35 35 | 65 2d 31 37 20 2b 20 30 |   -5.55|e-17 + 0|
|00000690| 69 0d 5c 65 6e 64 7b 76 | 65 72 62 61 74 69 6d 7d |i.\end{v|erbatim}|
|000006a0| 0d 0d 53 65 65 20 61 6c | 73 6f 3a 20 50 61 67 65 |..See al|so: Page|
|000006b0| 7e 5c 70 61 67 65 72 65 | 66 7b 66 6e 3a 65 69 67 |~\pagere|f{fn:eig|
|000006c0| 7d 20 5c 76 65 72 62 2b | 65 69 67 28 29 2b 2c 0d |} \verb+|eig()+,.|
|000006d0| 50 61 67 65 7e 5c 70 61 | 67 65 72 65 66 7b 66 6e |Page~\pa|geref{fn|
|000006e0| 3a 65 69 67 73 7d 20 5c | 76 65 72 62 2b 65 69 67 |:eigs} \|verb+eig|
|000006f0| 73 28 29 2b 2e 0d 0d 5c | 6e 65 77 70 61 67 65 0d |s()+...\|newpage.|
|00000700| 0d 0d 0d 0d             |                         |....    |        |
+--------+-------------------------+-------------------------+--------+--------+